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Abstract 

Elcctroosmotic pumping through uncharged hydrogcls can be achieved by embed- 
ding the polymer network with charged colloidal inclusions. Matos et al. (2006) re- 
cently used the concept to enhance the diffusion-limited flux of uncharged molecules 
across polyacrylamide hydrogel membraness for the purpose of improving the per- 
formance of biosensors. This paper seeks to link their reported macroscalc diagnos- 
tics to physicochemical characteristics of the composite microstructure. The exper- 
iments are characterized by a Debye screening length that is much larger than the 
radius of the silica nanoinclusions and the Brinkman screening length of the polymer 
skeleton. Accordingly, closed-form expressions for the incremental pore mobility are 
derived, and these are evaluated by comparison with numerically exact solutions 
of the full electrokinetic model. A mathematical model for the bulk electroosmot- 
ically enhanced tracer flux is proposed, which is combined with the electrokinetic 
model to ascertain the electroosmotic pumping velocity from measured flux en- 
hancements. Because the experiments are performed with a known current density, 
but unknown bulk conductivity and electric field strength, theoretical estimates of 
the bulk electrical conductivity are adopted. These account for nano-particle polar- 
ization, added countcrions, and non-specific adsorption. Theoretical predictions of 
the fiux enhancement, achieved without any fitting parameters, are within a fac- 
tor of two of the experiments. Alternatively, if the Brinkman screening length of 
the polymer skeleton is treated as a fitting parameter, then the best-fit values are 
bounded by the range 0.9-1.6 nm, depending on the inclusion size and volume frac- 
tion. Independent pressure-driven flow experiments reported in the literature for 
polyacrylamide gels without inclusions suggest 0.4 or 0.8 nm. The comparison can 
be improved by allowing for hindered ion migration, while uncertainties regarding 
the inclusion surface charge are demonstrated to have a negligible influence on the 
electroosmotic flow. Finally, and perhaps most importantly, anomalous variations 
in the flux enhancement with particle size and volume fraction can be rational- 
ized at present only by acknowledging that particle-particle and particle-polymer 
interactions increase the cfl'cctivc permeability of the hydrogel skeleton. This bears 
similarities to the increase in polymer free volume that accompanies the addition of 
silica nanoparticles to certain polymeric membranes. 
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1 Introduction 



Electric fields are widely used to probe the physicochemical characteristics 
of dispersed colloidal particles. Microelectrophoresis, capillary electrophoresis 
and gel electrophoresis, for example, are routinely used to characterize and 
sort DNA fragments, proteins, synthetic and biological macromolecules, and 
supramolecular aggregates. In all these applications, the electric field induces 
translation (electrophoresis) of charged particles, with the electrophoretic ve- 
locity reflecting a balance of electrical, hydrodynamic and steric interaction 
forces, which, in turn, are related to the particle size and charge, and properties 
of the dispersing medium. 



Matos et al. (2006) recently proposed a novel means of promoting the transfer 



of uncharged molecules across uncharged polyacrylamide hydrogels. By em- 
bedding the polymer skeleton with immobilized silica nanoparticles, they en- 
dowed the matrix with a fixed charge, and demonstrated that the application 
of an electric field enhanced the otherwise diffusion-limited flux of uncharged 
tracer across the membrane. Despite the possibility of promoting unwanted 
electrostatic interactions with charged molecules, the ability to control the ef- 
fective charge density by varying the concentration of immobilized inclusions 
offers several advantages over chemical and/or physicochemical modifications 
of the polymer architecture alone. Specific interactions may also be imparted 
by decorating the inclusions with functional moieties, such as in affinity chro- 
motography. 



As expected from studies involving electroosmotic microfiuidic pumping (Yao 



and Santiago 2003 Yao et al. , 2003), enhanced transport in immobilized cell 



culture (Chang et al. , 1995) and electroosmotic dewatering (Yukawa et al 



1976), the origin of the electric-field induced flux enhancement in Matos et 



al.'s experiments was qualitatively attributed to electroosmotic flow originat- 
ing from the diffuse layers of charge that envelop the silica inclusions. More- 



over, the microstructure of the composites synthesized by Matos et al. has 
potential to serve as a model system with which to test electrokinetic theory 
that predicts, from first principles, the average electroosmotic pumping veloc- 
ity, electrical conductivity, and effective ion diffusion coefficients in this novel 
class of nano-structured composites (Hill, 2006cp ). Indeed, the purpose of this 
paper is to provide a quantitative theoretical interpretation of the available 
experiments and, in turn, link macroscale diagnostics to characteristics of the 
microstructure. 



The paper is set out as follows. First, a detailed description of Matos et al.'s 
experiments is provided in §2. This identifies the relevant parameters and 
highhghts the available diagnostics. Section 3 presents various aspects of the 
theory. In particular, §§3.1 describes the electrokinetic model; §§3.2 derives the 
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means of linking Matos et al.'s flux enhancement ratio to the electroosmotic 
pumping velocity; §§3.3 addresses the bulk electrical conductivity of the com- 
posite membranes; and §§3.4 evaluates — by comparison with numerical solu- 
tions of the full electrokinetic model — a convenient closed-form approximation 
for the incremental pore mobility. The results are presented in §4. Here, lit- 
erature values of the permeability of gels without inclusions — obtained from 
pressure-driven flows — are compared to values inferred from the electrokinetic 
interpretation of the experiments. Certain anomalies are discussed in §5, and 
the paper concludes with a summary (§6). 



2 Experimental 



The hydrogel nanocomposite membranes synthesized in Matos et al.'s experi- 
ments comprise a polyacrylamide hydrogel matrix with silica nanosphere inclu- 
sions. The hydrogels were synthesized with a polymer mass fraction c = 0.15 
(mass of polymer to mass of water). After polymerization and cross-linking, 
however, the matrix was reported to swell, yielding a flnal polymer mass frac- 
tion c ~ 0.13. The microstructure of the membranes is characterized by the 
radius of the inclusions (a ~ 3.5 or 7.5 nm) and the inclusion volume fraction 
(0 ^ 0.02 or 0.04). 

A schematic illustration of Matos et al.'s apparatus is presented in flgure 1. 
The membrane comprises 29 cylindrical sub-membranes, each with a diameter 
2.68 mm (giving a total cross-sectional area A 1.64 cm^) and thickness L ^ 
1.68 mm. Matos et al. measured the diffusion coefficient of the tracer (amino- 
methylcoumarin) in gel membranes with and without inclusions, establishing 
a value D ^ 3.3 x 10~^° m^s~^. They did not discern the 0(0) -C 1 influence 
of the impenetrable inclusions on the effective diffusion coefficient. 

The experiments were initiated by establishing a quasi-steady diffusive flux of 
uncharged tracer across each membrane before applying a potential difference 
to the electrodes. The electric fleld was applied after 30 min, which is short 
compared to the characteristic diffusion time scale L^/D ^ 140 min. However, 
it is evident from the experiments and, indeed, the exact solution of the tran- 



sient diffusion problem (Matos et al. , 2006), that the fastest transients decay 
in about 30 min. 

The gold electrodes, which present a circular surface area on each side of the 
membrane, have a diameter 1.9 cm and are separated by a distance 6 mm. In 
the experiments considered in this paper, the differential voltage was controlled 
to maintain a constant electrical current i ~ 3 mA. 

The experiments were performed with the source and the sink reservoirs ini- 
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Fig. 1. Schematic of the experimental apparatus used by Matos et al. (20061 to mea- 
sure the electric-field enhanced flux of a neutral tracer across polyacrylamide hydro- 
gel-silica nanocomposite membranes. The immobilized silica nanoinclusions endow 
the microstructure with a negative charge, which, with an electric field E = — di/i/dx 
between the electrodes, drives electroosmotic flow from the anode left) toward 
the cathode (— , right). The flow produces a convective enhancement of the oth- 
erwise purely diffusive flux of the neutral tracer, whose concentration is denoted 
n°°{x), from the source resevoir to the sink resevoir. See the text and Matos et al.'s 
paper for details. 

tially filled with either 1 mmol 1~^KN03 electrolyte or deionized (DI) water. 
Only the source reservoir contained tracer (with concentration 0.1 mmol 1^^) 
at the beginning of the experiments, and the electrolytes on each side of the 
membrane were recirculated at a rate Q ~ 8 cm^ min~^, with respective source 
and sink reservoir volumes Vi ~ 50 cm^ and V2 ~ 10 cm^. 



The residence time of fluid in the two gaps between the membrane and the 
electrode is St ^ 3.6 s, so the thickness of a diffusive boundary layer at the 
membrane surface is 61 ~ V 6tD ^ 35 jjm. Equating the diffusive flux across 
the membrane DAn°°/L to the flux across the boundary layer D6n°°/6l in- 
dicates that 6n'^/An°° = 61/ L ~ 10^^^, so the tracer concentration at the 
edges of the membrane is approximately equal to the bulk concentration in 
the respective reservoir. Note that the hold-up time of the tracer in the reser- 
voirs is less than Vi/Q ~ 6 min, which is short compared to the ~ V2L/{AD) 
characteristic time scale for the (diffusive) tracer flux to change the reservoir 
concentration by an 0(1) amount. It follows that the membrane and reservoir 
dynamics are quasi-steady. 

The tracer (amino-methylcoumarin) concentration was measured with a fluo- 
rescence spectrophotometer that continuously sampled fluid in the sink reser- 
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voir circuit. Matos et al. also measured the pH in both reservoirs, showing 
that the source and sink reservoirs, respectively, attain acidic and basic pH 
with the application of an electric field. 

Note that the tracer has an acid dissociation constant Ka = 10~^ mol 1~^, so 
the fraction of it in a protonated state, 1/(1 + i^'a/PsO"^]), is significant only 
when pH < 5. These conditions prevail in the source reservoir and, possibly, 
throughout the composite membrane, suggesting that electrophoresis of the 
tracer should contribute to the net electric-field-induced flux enhancement. 
From the diffusion coefficient D ^ 3.3 x 10^^° m^ s~^, the electrophoretic mo- 
bility is (to a first approximation) eD/{kT) ^ 1.1 (/im s~^)/(V cm~^). There- 
fore, with a representative electric field ~ 10 V cm~^, the electrophoretic ve- 
locity of the tracer is Ut ~ 10 /xm s~^, which corresponds to Pe = UtL/D ~ 55. 
This is significantly greater than the values ascertained from electroosmotic 
pumping, suggesting that the fiux of protonated tracer might be dominated by 
electrophoresis. However, the experimental evidence is that the fiux is dom- 



inated by electro-osmosis (Matos et al. , 2006). It follows that pH > 5 in a 



significant portion of the composite hydrogel membranes. 



Finally, from the tracer concentration time series, [Matos et al.| reported av- 
eraged fiux enhancement ratios, which they defined as the ratio of the tracer 
concentration in the sink reservoir after 300 min (with an electric field applied 
after the first 30 min) to the concentration from an experiment without an 
electric field. They appropriately termed this ratio an enhancement factor, 
whereas the analysis undertaken in this paper approximates it as a flux en- 
hancement factor, J/Jq. The ratio of the fiuxes, which undoubtedly exhibit 
weak transients, are therefore assumed to be constant in this work. It is evi- 
dent that individual experiments refiect significant, non-systematic statistical 
variations, so the analysis below relies entirely on five averaged values of the 
reported enhancement factor, each with various particle radius a and volume 
fraction 0, and electrolyte (1 mmol I'^KNOs or deionized water), but fixed 
membrane geometry and electrical current. 



3 Theory 



3.1 Electrokinetic model and electoroosmotic flow in uniform composites 



The electrokinetic transport model adopted in this work is described in consid- 
erable detail elsewhere (Hill 2006d). Briefiy, it comprises the non- linear Pois- 
son Boltzman equation for electrostatics; electrolyte conservation equations 
involving diffusion, electromigration and convection; and Brinkman's equa- 
tions (with an electrical body force) for solvent mass and momentum conser- 
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vation (Brinkman, 1947). A solution of the equations, with a single spherical 
charged, impenetrable colloid embedded in an unbounded, rigid polymer gel 
yields the far-field (r oo) decays of the electrostatic potential ip, ion con- 
centrations Tij and fiuid velocity u. Of prime importance in this paper are 
perturbations to equilibrium that arise from a uniform electric field E. In the 
far-field, these take the following asymptotic forms: 



i;^-E -r + {E ■ er)D^r-^ (1) 
nj^nf + {E-er)Cfr-^ (2) 
u -2C^r~^{E ■ e,)e, - C^r~'\E ■ ee)ee. (3) 



Note that and eg are radial and tangential unit vectors, with the polar axis 
being colinear with E. The scalar coefficients , and C^rj/i"^ are, 
respectively, the dipole strengths of the electrostatic potential, ion concentra- 
tions, and pressure perturbations induced by the forcing E. In general, these 
so-called asymptotic coefficients are complicated functions of the particle ra- 
dius and C-potential, a and (, the electrolyte viscosity and Debye screening 
length, 7] and k~^, and the Darcy permeability of the polymer skeleton i"^ 
(square of the Brinkman screening length) . 



The single-particle problem leads to averaged equations that quantify the bulk 
ion fluxes and fluid flow in a dilute (0[1 + (/ta)"^]^ <^ 1) composite (Hill 



2006d). For the statistically uniform composites addressed in this work, the 



pressure and bulk ion concentration gradients are zero. Accordingly, with a 
uniform electric fleld and unidirectional flow and electrolyte and tracer fluxes, 
the average fluid velocity is 



U = -4>{3/a^)EC^, (4) 
where <C 1 is the inclusion volume fraction, and the average ion fluxes are 



J. = n'^U + ZjC-^nfE + 0(3/a=^)2,e-^nf + 0(3/a=^)D,ECf , (5) 

where Zj = ±z {j = 1, 2) are the electrolyte ion valences, Dj are the ion dif- 
fusion coefficients, kT is the thermal energy, and e is the fundamental charge. 
In Eqn. (5), the first two terms are, respectively, the convective and electromi- 
grative fiuxes that prevail in the absence of inclusions; and the last two terms 
are (microscale) electromigrative and diffusive contributions to the bulk flux 
arising from disturbances of the inclusions. 
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3.2 Uncharged tracer flux and flux enhancement 



The flux of uncliarged tracer is 

J = n'^U -DB + (f){3/a^)DC^B, (6) 



wliere = (l/2)a^ and B = dn°°/ dx is tlie bulk concentration gradient (Hill 



2006c) 



Clearly, the electric-field-induced flow U and electrolyte ion fluxes Ji and J2 
are decoupled from the tracer, but the tracer flux is influenced by the elec- 
troosmotic flow. Therefore, the tracer concentration satisfies the conservation 
equation 

= -dJ/dx, (7) 

where 

J = n'^U - D[l- {3/2)(t)]dn'^/dx. (8) 



Under quasi-steady conditions {dn°° /dt ^ 0), which will be justified shortly, 
the flux may be written in terms of the tracer concentration at the edges 
(x = and L) of the membran^^ 

J = f/[n-(x = L) - = 0)e^^]/{l - e^^). (9) 

Here, the Peclet number is 

Pe = UL/D = (f)MEL/D (10) 

and 

M = U/{E^) = -{3/a^)C^[Ka, (e/{kT), kE] (11) 
is termed the incremental pore mobility. 

Scaling arguments that also consider transients in the tracer concentration in 
the (well mixed) source and sink indicate that a quasi-steady approximation 
is reasonable when V/{LA) ^ 1 and Pe ~ 1; otherwise, 1 ^ Pe <^ V/{LA) 
or LA/V -C Pe ^ 1. Here, V is the smallest of the source and sink volumes. 
For Matos et al.'s experiments, ~ 10 cm^, A ~ 1 cm^ and L ~ 0.1 cm, so 
V/ (LA) ~ 10 3> 1, as required when Pe ~ 1. 

With zero electric field (Pe = 0), the tracer flux is 

Jo = -D[n°°(x = L) -n°°(x = 0)]/L, (12) 



The tracer concentration is n = n°°{x = 0)e /■^ + (J/C/)(l — e /^). 
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so the flux enhancement J/ Jq may be written 

J/ Jo = [1 + e(l - eP^)]Pe/(eP^ - 1), (13) 
where e = n°°{x = 0)/[n°°{x = L) — n°°{x = 0)]. 

If Jo > with e ~ — 1, which corresponds to diffusion from the source to the 
sink, with the tracer concentration in the sink much less than in the source, 
it may be reasonable to write 

J/Jo^Pe/(l-e-Pe). (14) 

Clearly, the flux is enhanced by an electric field that produces an electroosmotic 
flow directed down the concentration gradient (Pe > 0). 



3.3 Polarization, added counterions and non-specific adsorption 



The ion mobilities in the theory above are unhindered by the polymer skele- 
ton when they are calculated from widely available limiting conductivities. 
Therefore, the bulk electrical conductivity of the gel in the absence of in- 
clusions takes the same value as the electrolyte in the absence of polymer. 
Furthermore, the experimentally determined mobilities depend on a theoreti- 
cal prediction of the bulk electrical conductivity of the composite, taking into 
account concentration and electrical polarization, added counterions, and non- 
specific adsorption. Accordingly, the effective conductivity may be written 

= K^[1 + 0(Ap + A,, + A„,J], (15) 

where the polarization increment Ap is obtained from the ion fluxes above, 
giving 

Ap = i3/a')D^ey{kT) ^ zp^nf/K^ + (3/a=^)e ^ z^Dfif/K^. (16) 

i=l,2 i=l,2 

This reflects the electric-field-induced electrostatic potential and concentration 
dipole moments [D^ and D^) for a single immobilized inclusion (Hill 



2006d 



In contrast, the added counterion and non-specific adsorption increments do 
not depend on perturbations induced by the electric field; they can therefore 



be calculated from knowledge of only the equilibrium base state (Saville, 1983 



Hill et al. , 2003), as parameterized by (e/{kT) and Ka. 



The total conductivity increment Ap + Aac + Ansa is plotted as a function of 
Ka for various (^-potentials in figure 2 (solid lines). Recall, the added counte- 
rion and non-specific adsorption contribution Aac + Ansa is an 0(0) contribu- 
tion to the equilibrium base state, whereas the polarization contribution Ap 
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Fig. 2. The conductivity increment (3/2 has been added to permit a logarithmically 
scaled ordinate) A'^ = Ap + A^c + A„sa (solid lines) as a function of the reciprocal 
double layer thickness na for various (scaled) C-potentials —C,e/{kT) = 1, 2, 3, 
6: aqueous KNO3 at T = 25°C; a = 10 nm; 1 = 1 nm. Dashed lines are the contri- 
butions from electrostatic and concentration polarization, A^; these are practically 
independent of the Brinkman screening length I. The horizontal line at 3/2 denotes 
where A'^ = 0. 



(dashed lines) is an 0(0) contribution to the perturbed equilibrium. At low 
ionic strengths (small /ta), Ap is smaller than for dispersions where the parti- 
cles undergo electrophoresis in a Newtonian electrolyte (Hill, 2006d). There- 
fore, when the particles are immobilized in a polymer network, the added 
counterion and non-specific adsorption contribution plays a significant role. 



Note that the calculations presented below for the experiments with DI water 
as the electrolyte assume a C-potential that is the same as in 1 mmol 1~^ 
KNO3 (C = — 70 mV); the calculations also neglect the presence of dissolved 
CO2. The later can be justified from the bulk conductivity being dominated by 
the added counterions (H3O"'"), while the former is speculative and principally 
motivated by the unreasonably high ^-potentials that result when a constant 
surface charge is assumed. Nevertheless, when 0(Aac -|- A„sa) ^ 1, which is 
the case for DI water and 1 mmol 1~^ KNO3, is influenced by C, in the 
same manner as M. Therefore, while the value of C, chosen to estimate 
influences the experimentally ascertained pore mobility M [see Eqn.(19)], it 
affects the theoretical prediction of M in the same manner, thereby yielding 
a similar best-fit value of £. 
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3-4 Asymptotic formulas for the pore mobility 



In addition to the numerically exact solutions of the electrokinetic model, let 
us recall an approximate analytical solution valid for dilute composites where 
the ^-potential is low (|C| < kT/e) and the double-layer thickness and 
Brinkman screening length i are both small compared to the particle radius 
a: 

This formula confirms inferences drawn from numerical solutions of the full 



electrokinetic model (Hill, 2006d), and, indeed, Matos et al.'s experiments, 
namely that the electroosmotic flow is inversely proportional to the particle 
radius. However, Matos et al.'s experiments are characterized hj na < 1 with 
£ ~ a. As shown below, such conditions produce an electroosmotic flow that 
scales with the reciprocal square of the particle radius at constant ^-potential 
and electric field strength. 

To obtain a convenient formula for small na, note that the hydrodynamic 
force on an immobilized charged colloid vanishes with respect to the electrical 



force when /ta -C 1 and \(\ < kT/e (Hill, 2006a). Further, when Ka — >■ and 
Id < kT/e, the electrical force approaches the bare Coulombe force f^ — 



AneseoaCE. Therefore, since = ATT{ri/f)C^E ( [HiUl|2006c[ ), the incremental 
pore mobility becomes 

f//(E0) = -3e,eoC(^/a)V^ (f^a < 1, |C| < kT/e). (18) 

It is noteworthy that the dependence of the pore mobility on particle size is 
qualitatively different at high and low na. This is due to the relationship be- 
tween the surface charge density and potential. When na is large, the surface 
charge density is ~ (es€oK,, so the average mobile charge density in the hydrogel 
scales as (eseoH(f)/a. However, when na is small, the surface charge density is 
~ (egto/a, so the mobile charge density scales as (eseo(p/a'^- Accordingly, bal- 
ancing the electrical and Darcy drag forces gives — -E47r^eseoa0/[(4/3)7ra^] = 
{ri/i'^)U, so the incremental pore mobility is M = U/{E(f)) = —?)eseoC{(-/a)'^/rj, 
which agrees with Eqn. (18). 



Matos et al. (2006) identified the electroosmotic flow in their experiments as 
being inversely related to the surface area density (p/a. However, their data 
was insufficient to ascertain a specific power-law dependence. While the theory 
clearly identifies the mobility as having a reciprocal square dependence with 
particle size at constant ^-potential when <^ 1, and a reciprocal dependence 
when Ka ^ 1, the theoretical interpretation of the experiments is further com- 
plicated by the conductivity of the composite depending on the particle size 
and volume fraction. In fact, when na <^ 1, / = 1 + O[0(Ka)~^], so 
U = (pME has a weaker dependence on particle size and volume fraction 
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10'^ lo' 10° lo' 10^ lo' " 10° lO' 10^ lo' 



scaled reciprocal double layer thickness, Ka scaled reciprocal double layer thickness, Ka 

Fig. 3. The incremental pore mobility M = U/{E(j)) as a function of the scaled 
reciprocal double layer thickness na for various (scaled) ("-potentials —C,e/{kT) = 1, 
2, 3, 6: aqueous KNO3 at T = 25°C; a = 3.5 (left panel) and 100 nm (right panel); 
£ = 1.0 nm. Solid lines are numerically exact solutions of the full elcctrokinctic 
model, and the dashed lines are approximate analytical solutions [Eqns. (17) (large 
Ka) and (18) (small Ka)]. 

than expected if the effective conductivity were assumed constant. More im- 
portantly, and, perhaps, surprisingly, it will be demonstrated that {{) at very 
low bulk ionic strengths the electroosmotic pumping velocity U at constant 
current density I — i/ A is independent of the particle size, charge and vol- 
ume fraction; and (ii) the effective Darcy permeability of the polymer skeleton 
evidently increases with decreasing average interparticle surface proximity. 

Representative numerical calculations of the pore mobihty are compared with 
the asymptotic formulas in figure 3 for inclusions with two representative radii 
(a — 3.5 and 100 nm) dispersed in an hydrogel with Brinkman screening length 
£ = 1 nm. In regions of the parameter space where the full model agrees with 
Eqn. (18) (small ko) and (17) (large kg with i ^ a), mobilities for a different 
values of i can be conveniently obtained by multiplying the ordinate by £^ 
(measured in nm^). 



4 Results 

The principal methodology adopted to interpret Matos et al.'s experiments is 
summarized as follows. First, the measured flux enhancements and Eqn. (14) 

are used to determine Pe. Then, from knowledge of the electrical current {i = 
3 mA), membrane geometry {A = 1.64 cm^ and L = 1.68 mm), particle 
volume fractions (0 = 0.02 and 0.04) and radii (a = 3.5 and 7.5 nm), and a 
theoretical calculation of the effective conductivity K^, the incremental pore 
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mobilities are 

M = FeD/{(j)EL) = VeDK''A/{(t)iL). (19) 

With these 'experimentally determined' mobilities, the electrokinetic model 
provides the effective hydraulic permeabilities £^ of the hydrogel in each com- 
posite. In this manner, the theoretical interpretation, which necessitates theo- 
retical predictions of and M, is evaluated by the consistency of the resulting 
Brinkman screening lengths £. 

The electrokinetic transport calculations presented in table 1 were performed 
with unhindered ion mobilites. However, the conductivity increment and in- 
cremental pore mobility are practically independent of the ion mobilities, so 
it is reasonable to account for hindered ion mobilities by multiplying the un- 
hindered bulk conductivity by a single hinderance factor. This increases the 
calculated electric field strength and, therefore, decreases the incremental pore 
mobility inferred from the experimentally measured flux enhancement. In turn, 
smaller best-fit values of the Brinkman screening length emerge. 

The table summarizes theoretical interpretations of the five experimental con- 



ditions for which Matos et al.l reported statistically averaged flux enhance- 



ments. Cases A-C and D-E, respectively, identify experiments with 1 mmol 1~^ 
KNO3 and DI water. Note that DI water is approximated here as an elec- 
trolyte comprised of H3O"'" and OH^ ions, each with a concentration n°° = 
10~^ mol 1~^, and the (^-potentials of the silica nanoparticles are taken to be 



—70 mV in 1 mmol 1 KNO3 and DI water; recall, Matos et al. reported 
C-potentials only for 1 mmol 1~^ KNO3. 

The upper section of the table lists quantities that are readily available from 
Matos et al.'s paper. In addition to the flux enhancement, J / Jo, is a length 
scale a(0~^/^ — 2) that characterizes the average distance between the in- 
clusion surfaces. The middle section presents theoretical calculations of the 
polarization and added counterion and non-specific adsorption contributions 



to the conductivity increment (Saville, 1983). With the inclusion volume frac- 



tions and bulk electrolyte conductivity these provide the bulk conduc- 
tivity of the composite K'^. In turn, knowledge of the current density X = 
i/Ak. 1.83 mA cm~^ (in the composite) provides the electric field strength 
E = X/K^. Note that the conductivities of the composites saturated with 
DI water are comparable to those with 1 mmol 1~^ KNO3, even though the 
bulk ionic strength and, hence, conductivity of the electrolyte is four orders 
of magnitude lower. Clearly, particle polarization and the added counterions 
play a very important role. 

The lower section of the table presents calculations of the 'experimentally 
determined' Peclet number Pe = UL/D and incremental pore mobility M = 
U / {E(j)). Note that Pe may be considered a dimensionless electroosmotic pump- 
ing velocity here, since the membrane thickness L and tracer diffusivity D are 
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same in all cases A-E. The table presents theoretical predictions (from numer- 
ical solutions of the full electrokinetic model) of the incremental pore mobility 
for composites with a Brinkman screening length i = 1 nm. These have been 
used to calculate best-fit values of i based on knowledge that M is inversely 
proportional to i"^ at the prevailing values of < 1 (see figure 3). 

The scenario that emerges from these calculations is, perhaps, more complex 
than suggested by Matos et al.'s qualitative interpretation of the data. For 
example, comparing the Peclet numbers for cases A and B, where the particle 
size is doubled at a fixed inclusion volume fraction, suggests that U scales as 
a~^. While this is consistent with theoretical expectations if £ is a constant, 
doubling the particle size decreases the mobility by almost an order of mag- 
nitude. The anomaly can be reconciled by allowing i to vary, but it is not 
clear that the resulting best-fit values of i are acceptable. A similar discussion 
applies to the comparison of cases D and E (with Dl water). Next, comparing 
cases B and C, where the inclusion volume fraction is doubled with a fixed par- 
ticle size, the Peclet numbers indicate that U more than quadruples, whereas 
the theory requires U to double. Furthermore, the mobility triples, while the 
theory demands a constant mobility if i is constant. 

Figure 4 compares the Brinkman screening lengths ascertained from inde- 
pendent pressure-driven flows in thin polyacrylamide gel membranes without 



nano-inclusions. The two sets of data reported by White (1960) and Tokita 



and Tanaka (1991) are clearly very different, and it is unfortunate, perhaps, 
that the later authors did not reconcile these discrepancies. Nevertheless, in 
favor of their results is the consistency of their data with scaling theory. Ac- 
cordingly, a semi-empirical fit to their data, 

i ^ 0.182c"°-^35 ~ (nm), (20) 



gives i ~ 0.8 nm at the polymer concentration c = 0.13 reported by Matos 



et al. ( 2006[ ). This is clearly smaller than all the values of i (spanning the range 



0.93-1.62 nm) ascertained from the theoretical interpretation of Matos et al.'s 
data presented in table 1. 

Note, with a fixed value of i = 0.8 nm, table 2 presents theoretical calcula- 
tions with parameters corresponding to cases A-E in table 1. Again, all the 
calculations (leading to the incremental pore mobility and electrical condic- 
tivity) have been performed with unhindered electrolyte ion mobilities. The 
results in the upper and lower sections of the table are, respectively, with a 
1 mmol 1~^ KNO3 electrolyte (with K"*" counterions) and DI water (with HsO"*" 
counterions) . 

The calculations with a KNO3 electrolyte indicate that the mobilities, Peclet 
numbers, and flux enhancements exhibit similar qualitative trends as the ex- 
periments. Howevever, the quantitative trends are poor. Note that the con- 
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Fig. 4. The Brinkman screening length i (nm) of polyacrylamide hydrogels, as- 
certained from pressure-driven flow experiments in thin membranes, as a function 



of the polymer (mass) fraction c: White (19601 (squares) and Tokita and Tanaka 



(1991 1 (circles). Curves are empirical fits to the data: i = 0.182c '^■^^^ (nm) (solid); 



^ = 0.735 exp (-4.36c) (nm) (dashed). 

ductivity increments are practically the same as those presented in table 1 
with i = 1 nm. This reflects the extremely weak role of microscale convec- 
tion in polarizing the inclusions, and, of course, that added counterions and 
non-speciflc adsorption are an 0(0) contribution to equilibrium. 

The calculations in the lower section of the table, with DI water as the elec- 
trolyte, demonstrate that the electroosmotic flow velocity, as indicated by Pe 
or J/Jq, is practically independent of the particle size and volume fraction 
at sufficiently low bulk ionic strengths. If the discrepancy between the elec- 
troosmotic flow velocity U = PeD/L ^ 37 nm s~^ and the experimentaly 
determined value U 74: nm s~^ for DI water (when a(0~^/^ — 2) 11.8) 
can be attributed to the neglect of hindered ion migration, which leads to an 
incorrect estimate of the bulk conductivity K°°, then the effective electrolyte 
ion hinderance factor is suggested to be about 37/74 ^ 0.50. This is in rea- 
sonable agreement with independent measurements of the hindered diffusion 
coefficients of KCl and D2O (and two larger molecules) in polyacrylamide hy- 
drogels. In particular, at the bulk polymer concentration used in Matos et al.'s 



experiments (c ~ 0.13), White and Dorion (1961) measured diffusion coeffi- 



cients of KCl and D2O that are factors of 1.2/1.8 ^ 0.67 and 1.2/2.2 ^ 0.55 
smaller than their respective unhindered values. 

Accordingly, if it is assumed that the electrical conductivity of the membranes 
is indeed reduced by a factor of 0.50, and, furthermore, the mobilities of all 
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electrolyte ions are reduced by the same factor, then the incremental pore 
mobilities M inferred from Matos et al.'s experiments should be a factor 0.50 
smaller than the values listed in table 1. Moreover, the best-fit Brinkman 
screening lengths ^ in the table should be reduced by a factor v^O.50 ~ 0.71. 
These values are denoted ^h = V0.50£ in the lower section of table 1. It is 
noteworthy that their range 0.66-1.14 nm is much more representative of the 



values expected from the pressure-driven flows of White (1960) and Tokita 



and Tanaka (1991). 



Knowledge of the electrical conductivity of the electrolyte-saturated hydrogel 
would provide an independent quantitative measure of the influence that the 
polymer has on electrolyte ion mobilities while under the influence of an elec- 
tric field. Unfortunately, the electrical conductiviy of the hydrogels in Matos et 



al.'s experiments is unknown. However, data in Matos's thesis (fig. 4.29 Matos 



2006) indicates that the voltage drop across the entire cell {i.e., between the 
electrodes) when i = 3 mA is about 12 V when gels without inclusions are 
saturated with 1 mmol 1~^KN03. If this voltage is assumed to reflect only the 
combined electrical resistances of the membrane and electrolyte sandwiched 
between the electrodes, then the conductivity of the membrane is 

= K'^/[K°^AA/{tL) - {L'/L){A/A')]. (21) 

Here, A ^ 12 V is the total potential drop, K'^ ^ 14.5 mS m~^ is the 
conductivity of 1 mmol 1~^KN03 electrolyte (from knowledge of the limit- 
ing conductances of and NO3 ), L' = 6 — 1.68 = 4.32 mm is twice the 
distance between the membrane and one electrode, and A' f» 2.84 cm^ is 
the clear cross-sectional area between the membrane and electrodes. It fol- 
lows that K'' ^ ir°°/(5.66 - 1.56) ^ 0.24ir°°. While this indicates that the 
polymer significantly hinders ion migration, the calculation is likely to over- 
estimate the hindrance, because it neglects the (unknown) potential drop due 
to electrochemical reactions. 

Nevertheless, with a hinderance factor of 0.24 adoped to provide a lower- 
bound on the membrane electrical conductivity, the corresponding Brinkman 
screening lengths, denoted ih = \/0.24£ in the lower section of table 1, are 
a factor VO.24 ^ 0.49 smaller than the best-fit values ascertained with un- 
hindered ion mobilities. The range ih = 0.53-0.79 nm is representative of the 



values expected from the pressure- driven flows of White (1960) and Tokita 



and Tanaka (1991), but, as highlighted below, such a favorable comparision 



may be misleading. 

The analysis summarized in table 1 has led to an upper range i ~ 0.9-1.60 nm 
ascertained with unhindered ion mobilities, and a lower range ih ~ 0.53- 
0.79 nm ascertained using a lower-bound estimate of hindered ion mobilities. 
Note that fluctuations within each range correspond to unreasonably large 
fluctuations in the bulk polymer density, and such fluctuations cannot be rec- 
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onciled by osmotic swelling. It therefore remains to elucidate these significant 
variations in gel permeability. 



5 Discussion 



One clue to resolving the apparent shortcomings is to recognize that two dis- 
tinctly different values of the best-fit Brinkman screening lengths emerge from 
the experiments. Considering, for example, the values ascertained with unhin- 
dered ion mobilities, cases B and E yield values of I close to 1 nm, whereas 
cases A, C and D yield values close to 1.5 nm. These differences do not cor- 
relate with the electrolyte, particle size or volume fraction alone. However, 
there is strong correlation with a length scale that characterizes the average 
distance between closest neighbors in a uniform random dispersion. Letting 
the characteristic center-to-center distance be a0~^/^, the distance between 
the surfaces is a(0~^/^ — 2). This measure is listed in the last column of the 
upper section of table 1, and its relationship to the effective permeability is 
plotted in figure 5. 

Recall, the best-fit values of £, ascertained without hindered ion mobilities in- 
fiuencing the membrane electrical conductivity, provide an upper-bound on the 
effective gel permeability (open symbols). Furthermore, the lower-bound on 
the electrical conductivity, established above from the voltage-current relation- 
ship for the cell, provides a lower-bound-estimate of the gel permeability (filled 
symbols). As demonstrated in the figure 5, these upper- and lower-bounds ap- 
proach values expected from independent pressure-driven fiow experiments 



in homogeneous (i.e., without inclusions) polyacrylamide gels (Tokita and 



Tanaka, 1991 White, 1960). Moreover, the approach occurs as the separation 



a(0 ^'^ — 2) becomes large compared to 

This quantitative analysis supports Matos et al.'s suggestion that particle 
interactions play a necessary role in endowing the composites with a significant 
electroosmotic pumping capacity. Nevertheless, the physical origin of such 
interactions and, indeed, their quantitative infiuence on the electroosmotic 
pumping capacity and bulk conductivity remain unknown. 



Matos et al.| suggested that electroosmotic fiow creates pathways that some- 
how link particles in close proximity, as expected in a percolating network. 
However, their neutron scattering data did not reflect fractal-like correlations. 
Another possibility is that polymer structure {e.g., partial chain orientation 
and/or cross-linking density) is altered near the particle surfaces, producing 
shells of hydrogel that are much more electrically conductive {i.e., with higher 
ion permeabilities) and hydrodynamically permeable than in the bulk. Swelling 
after the synthesis could also generate shells of void space. 
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interparticle separation, a(i|) ^ -2) (nm) 



Fig. 5. The effective Brinkman screening lengtli as a function of the average 
distance between the surfaces of nearest neighbor inclusions in random disper- 
sions a((/)~^/^ — 2): upper-bound estimates based on theoretical interpretation of 
experiments with unhindered ion mobilities (open symbols) ; lower-bound estimates 
achieved with a lower-bound estimate of hydrogel electrical conductivity (filled sym- 
bols). See text for details. Circles and squares denote, respectively, cases A-C with 
1 mmol 1~^KN03 electrolyte, and cases D and E with DI water (see table 1). The 
solid lines are empirical best-fit (quadratic) decays: = £[l-|-47.3/{((/)~^/'^— 2)a/^}^] 
with I = 0.83 nm (upper curve); t = i{\^\Ml{{(\)-^l'^ -2)al 1^] with I = 0.42 nm 
(lower curve). Dashed lines indicate the values of I expected from pressure-driven 
flows in polyacrylamide gels with a polymer concentration c ~ 0.13: I ~ 0.8 nm 
( [Tokita and Tanaka[ |1991[ ) (upper hne); I 0.4 nm ( |White[ |1960| ) (lower line). 



It is tempting to draw upon Matos et al.'s data with a = 3.5 nm and (p = 0.005, 



0.01 and 0.2 (as seen in fig. 8 of Matos et al. , 2006) to help elucidate the 



role of particle volume fraction and size. However, it is evident from these 
single experiments, and the averaged data (used exclusively above), that there 
are significant statistical fluctuations from one experiment to another. For 
example, the flux enhancements from the experiments presented in figures 1 
and 8 of Matos et al.'s paper with DI water (with a = 3.5 nm and = 
0.02) are, respectively, ~ 1.58 and 1.63, whereas the flux enhancement from 
repeated measurements under the same conditions was reported as 1.44 ±0.01 
(as adopted in table 1 above). It is also difficult to draw firm conclusions from 
the experiments with DI water, since the ^-potentials are unknown, further 
hampering theoretical estimates of the bulk conductivity and incremental pore 
mobility. 



Finally, note that silica nanoinclusions have been demonstrated to enhance the 



permeability of glassy polymeric matrices to gases (Merkel et al. , 2002). This 
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unexpected behavior is attributed to an overall increase in the free volume of 
the polymer. Moreover, recent theoretical studies show that the increase in 
free volume is localized to the particle-polymer interface (Hill, 2006b). The 
interpretation of Matos et al.'s experiments summarized in figure 5 suggests 
that a similar influence prevails when silica nanoparticles are immobilized in 
swollen water-saturated polymer networks. 



6 Summary 



The full electrokinetic transport model proposed by Hill (2006d) was used 



to interpret recent experiments where Matos et al. (2006) immobilized silica 



nanospheres in an uncharged polayacrylamide hydrogel matrices and applied 
an electric field to electroosmotically enhance the otherwise a purely diffusive 
flux of an uncharged tracer across the composite membrane. 

Convenient simplifications of the full model, leading to simple closed-form for- 
mulas, were derived for the limit where the bulk ionic strength is low and the 
inclusion radius is small; these complement earlier analytical solutions for sit- 
uations where the Debye and Brinkman screening lengths are small compared 
to the particle radius. The analytical results were demonstrated to yield accu- 
rate approximations of the incremental pore mobility from numerically exact 
solutions of the full electrokinetic model. Figure 3 identifies practical bounds 
on the parameter space where these convenient approximations are accurate. 

A simple (quasi-steady) convective-diffusion model for the tracer flux was 
derived to quantitatively interpret the fiux-enhancement diagnostic reported 



by Matos et al. The general solution depends on the (changing) concentration 
of tracer in the source and sink reservoirs, but it was assumed here that the 
concentration in the sink was small enough to permit an approximation that 
relates the scaled electroosmotic flow velocity (Pe = UL/D > 0) directly to 
the flux enhancement J/Jq>1. 

Several approaches were taken to compare the theory and experiments. One 
involved selecting the Brinkman screening length of the polymer skeleton as 
a fitting parameter and using the experiments and theory to establish values 
that yield 'correct' flux enhancements. With unhindered electrolyte migration, 
the effective Brinkman screening length of the hydrogels was found to span the 
range 0.93-1.63 nm, depending, foremost, on the average particle separation, 
as measured by a(0^^/^ — 2). Evidently, as the particle separation diverges, 
the effective Brinkman screening length approaches a value (~ 0.8 nm) that 



is expected from independent pressure-driven flow experiments by Tokita and 



Tanaka (1991). On the other hand, allowing for hindered electrolyte ion mi- 



gration, a lower-bound-estimate of the efffective Brinkman screening length 
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gives a range 0.46-0.79 nm. Similarly, as the particle separation diverges, the 
effective Brinkman screening length approaches a value (~ 0.4 nm) that is ex- 



pected from independent pressure-driven flow experiments by White (1960). 



Several means of reconciling the apparent differences between theory and ex- 
periment were proposed, all of which require more extensive and complete 
experimental investigation to draw firm conclusions. Firstly, this analysis sup- 
ports Matos et al.'s suggestion that particle interactions, which are neglected 
entirely in the present theory (due to the small inclusion volume fraction), 
significantly increase the electroosmotic pumping capacity of the composites. 
A model to quantify the influence of such interactions has not been pursued, 
and it is not clear that the enhanced permeability arises from particle-particle 
or particle-polymer interactions, or, perhaps, a combination of both. 

While the electroosmotic pumping capacity of a composite is independent 
of the bulk electrical conductivity, when the electric field driving the flow is 
known, in practice it is the electrical current that is known accurately. This 
necessitates an accurate experimental determination of the bulk conductivity 
to complement a satisfactory theoretical interpretation of experiments, and 
thus permit an accurate prediction of the pumping capacity as a function of 
current density. 

Clearly, it would be advantageous if future experimental efforts reported the 
conductivity of hydrogels with and without inclusions, preferentially when uni- 
formly saturated with well-characterized electrolytes, including DI water. Such 
diagnostics should be performed using a high-frequency electric field while the 
membrane is in its steady or quasi-steady state, i.e., when subjected to a 
steady electric field. Note, however, that care must be taken to separate the 



influences of the bulk and near-electrode regions ( HoUingsworth and Saville 



2003) 
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